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Abstract. The Bethe-Salpeter equation in a strong magnetic field is studied for 
positronium atom in an ultra-relativistic regime, and a (hypercritical) value for the 
magnetic field is determined, which provides the full compensation of the positronium 
rest mass by the binding energy in the maximum symmetry state. The compensation 
becomes possible owing to the falling to the center phenomenon. The relativistic 
form in two-dimensional Minkowsky space is derived for the four-dimensional Bethe- 
Salpeter equation in the limit of an infinitely strong magnetic field, and used for finding 
the above hypercritical value. Once the positronium rest mass is compensated by 
the mass defect the energy barrier separating the electron-positron system from the 
vacuum disappears. We thus describe the structure of the vacuum in terms of strongly 
localized states of tightly mutually bound (or confined) pairs. Their delocalization for 
still higher magnetic field, capable of screening its further growth, is discussed. 


1. Introduction 

In the present paper we are considering, in the framework of quantum electrodynamics, 
the phenomenon of falling to the center in a system of two charged particles caused by the 
ultraviolet singularity 1/x 2 of the photon propagator - which mediates the interaction 
between them - on the light-cone, x 2 = — x 2 ~ 0. Here Xq is the time, and x is 

the space coordinate. This phenomenon occurs in a number of problems. In some of 
cases listed below the matter comes to the one-dimensional Schrodinger equation with 
the singular attractive potential U(r ) = —/3/r 2 , 0 < r < oo, the falling to the center 
depending on the value of the constant [3 that fixes the strength of the coupling [TJ. 

Long ago it was established that the falling to the center is inherent in the Bethe- 
Salpeter equation for electron-positron system treated in a fully relativistic way and 
taken in the ladder approximation, provided that the fine structure constant a exceeds 
a critical value a cr , (1/137) -C a CI < 1. We refer to the so-called Goldstein solution 
[2j of the Bessel differential equation in the variable s = \/—x 2 , the space-like interval 
between the particles, to which the problem is reduced in the ultra-relativistic limit. 
Falling to the center implies that the singular attraction is so strong that the energy 
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spectrum becomes unlimited from below, and hence - once the singular problem is 
duely defined - the binding energy may completely cancel the rest mass, so that the 
gap between positive and negative continua disappears. This reminds to some extent 
the situation in the problem of an electron in an external electric field, occupying the 
whole space, where the tunnelling through the gap leads to the vacuum instability with 
respect to e + e _ — pair creation known as Schwinger effect, although the latter does not 
have a threshold character. The value a = a CI may delimit the range of self-consistency 
of the theory known as quantum electrodynamics, at least in its customary form. 

Another well-known theater 0) 0) 0 where the falling to the center acts is the 
point-like nucleus with its charge Z exceeding the critical value Z cr = a = 137, into 
whose Coulomb held a relativistic electron described by the Dirac equation is placed. 
After this problem is properly regularized, the electron energy level reaches the border 
e = —m of the lower continuum as the charge Z grows, which gives rise to electron- 
positron pair production, the free positrons leaving the atom. This mechanism leads to 
diminishing the charge and restores it at the critical value that cannot thus be exceeded. 

If an electron in an atom is treated within the Schrodinger equation, or a pair 
of particles is addressed in a semi-relativistic way with the aid of the Bethe-Salpeter 
equation, wherein the famous Salpeter’s equal-time Ansatz is made (resulting in the 
disregard of the retardation effects), the falling to the center does not take place. The 
matter is that the singularity in the point r = 0, r = of the Coulomb potential 
U(r ) = —aZ/r originating in this case from the light-cone singularity of the photon 
propagator, mentioned above, is not sufficiently strong, and the energy level remains 
shallow. On the contrary, in the above two cases the falling to the center may 
be attributed to relativistic enhancement of the Coulomb force - the second-order 
differential equation to which the spherically symmetric Dirac equation may be reduced 
by excluding the second spinor component, contains the sufficiently singular attractive 
term —( aZ/r ) 2 in the effective potential [B]. 

The situation changes drastically when a strong external magnetic held B is 
applied to the system. Already the Schrodinger equation for a particle situated in 
the combination of the Coulomb and strong magnetic fields possesses the falling to the 
center caused by the singularity l/\z\ in the differential Schrodinger equation defined 
on the whole axis — oo < z < oo. Here z is the electron coordinate along the axis 
parallel to the magnetic held. The energy spectrum of this equation contains large 
negative values tending to minus infinity jH] as the magnetic held B —> oo. The reason 
lies in the dimensional reduction: for strong magnetic held the electron is restricted to 
the lowest Landau level, consequently its wave function obeys a two-dimensional (one 
space and one time coordinates) differential equation. The reduction of the number 
of degrees of freedom causes the effective strengthening of the attraction, like in the 
quark confinement problem. Analogous strengthening takes place when the electron- 
positron system is described by the Bethe-Salpeter equation studied using the semi- 
relativistic Ansatz that implies the non-relativistic character of the relative motion of 
the constituent particles 0-mi Again, the energy spectrum of the corresponding two- 
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dimensional equation is unlimited from below as the magnetic field grows. The total rest 
mass of the ground state is compensated when the magnetic held reaches the value [9], 


whose order of magnitude is determined by the large factor exp (const /a). As pointed 


in {9] those results, however, are not reliable, since they depend on the unrighteous 
extrapolation of non-relativistic procedure to the essentially relativistic region of large 
mass defect. 

The relativistic treatment of an electron in a combination of external Coulomb and 
magnetic held was given in 03- The result was that if the magnetic held is reasonably 
strong (~ 10 15 G) the critical value of the nucleus charge Z may correspond to stable 
elements with Z 90. 

In the present paper we consider the system of two charged relativistic particles 
- especially the electron and positron - in interaction with each other, when placed in 
a strong constant and homogeneous magnetic held B. To this end we use the Bethe- 
Salpeter equation in the ladder approximation without exploiting any non-relativistic 
assumption. We derive its limiting - when B —> oo - form, which comes out to be a 
Bethe-Salpeter equation in two-dimensional Minkowsky space-time, covariant under the 
corresponding Lorentz subgroup - the boost along the magnetic held. Stress, that the 
two-dimensionality holds only with respect to the degrees of freedom of charged particles, 
while the photons remain 4-dimensional in the sense that the singularity of the photon 
propagator is determined by the inverse d’Alambert operator in the 4-dimensional, 
and not two-dimensional Minkowsky space. (Otherwise it would be weaker). The 
term responsible for interaction with an arbitrary electric held directed along B is 
also included and does not lay obstacles to the dimensional reduction. Throughout 
the paper we set H — c — 1 and use the Heaviside-Lorentz units, where a = e 2 /47r. 

^ I Gaussian B | Heaviside—Lorentz ' 

We make sure that in the case under consideration the critical value of the coupling 
constant is zero, a cr = 0, i.e., the falling to the center is present already for its genuine 
value a = 1/137, in contrast to the no-magnetic-held case, where a CI > 1/137. If the 
magnetic held is large, but finite, the dimensional reduction holds everywhere except 
a small neighborhood of the singular point s = 0, wherein the mutual interaction 
between the particles dominates over their interaction with the magnetic held. The 
dimensionality of the space-time in this neighborhood remains to be 4, and its size is 
determined by the Larmour radius Lb = (eB)^ 1 ^ 2 that is zero in the limit B = oo. The 
latter supplies the singular problem with a regularizing length. The larger the magnetic 
held, the smaller the regularizing length, the deeper the level. We find the value of the 
magnetic held - we call it first hypercritical field - 



( 1 ) 


where Ce = 0.577 is the Euler constant, that provides disappearance of the center- 
of-mass energy of the electron-positron pair and of its center-of-mass momentum 
component along B. We refer to this situation as a collapse of positronium. 

In discussing the physical consequences of the falling to the center we appeal to the 
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approach recently developed by one of the present authors as applied to the Schrodinger 
equation with singular potential m and to the Dirac equation in supercritical Coulomb 
held |5J. Within this approach the singular center looks like a black hole. The solutions 
of the differential equation that oscillate near the singularity point are treated as free 
particles emitted and absorbed by the singularity. This treatment becomes natural 
after the differential equation is written as the generalized eigenvalue problem with 
respect to the coupling constant. Its solutions make a Hilbert space and are subject to 
orthonormality relations with a singular measure. This singularity makes it possible for 
the oscillating solutions to be normalized to 5-functions, as free particle wave-functions 
should be. The nontrivial, singular measure that appears in the definition of the scalar 
product of quantum states in the Hilbert space of quantum mechanics introduces the 
geometry of a black hole of non-gravitational origin and the idea of horizon. The 
deviation from the standard quantum theory manifests itself in this approach only 
when particles are so close to one another that the mutual Coulomb held they are 
subjected to falls beyond the range, where the standard theory may be referred to as 
firmly established [5]. 

Within this approach the regularizing length provided by the Larmour radius is 
dealt with not as a cut-off, but as a lower border of the normalization volume, the event 
horizon in a way. Although the result o is obtained following the concept of Refs. 
[13], .5J, it can be reproduced without essential alteration within the standard cut-off 
philosophy, too. 

The most intriguing question is what happens after the magnetic held exceeds 
the first hypercritical value 0 - The solution of the Bethe-Salpeter equation in two- 
dimensional space-time in the ultra-relativistic limit studied in the present paper 
corresponds to formation of special ’’confined” states in the kinematical domain called 
sector III in h?, . [5j. (Within the standard approach these would be bound states, 
although this is less adequate). As the corresponding overall energy and momentum 
of such e + e~-state is zero, it is not separated from the vacuum by an energy barrier. 
Besides, this state is the one of maximum symmetry in the coordinate and spin space. 
Hence, it may be thought of as relating to the vacuum, as well, and describing its 
structure. The confined particles cannot escape to infinite distance from one another, 
on the contrary the probability density of the confined state is concentrated near the 
point s = 0, behind the horizon - as distinct from the ordinary bound state. 

The situation is expected to change as the magnetic held goes on growing. At a 
certain stage - we reserve the name second hypercritical for the corresponding value 
of the magnetic held - deconhnement of the above strongly localized states may occur. 
The corresponding solutions to the Bethe-Salpeter equation are not yet strictly obtained, 
which makes us describe the deconhnement more hypothetically. After the level deepens 
further, the center-of-mass 2-momentum gets into sufficiently far space-like region, and 
solutions oscillating at large distances between electron and positron appear. Thus, the 
state delocalizes. The delocalized electron-positron pairs produced from the vacuum, 
each particle on a Larmour orbit, should screen the magnetic held and stop its growing 
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above the second hypercritical value, this value being the absolute maximum of the 
magnetic field admitted within quantum electrodynamics. Simultaneously, the space¬ 
like total momentum provides the lattice structure to the vacuum. 

This resembles the case of supercritical nucleus where there are states (that belong 
to sector IV in the nomenclature of Refs. [ 13 ], 0 admitting the leakage to infinity, which 
provides the mechanism for reducing the charge of the nucleus below the critical value. 

No sooner than the delocalized states are found in our present problem one may 
definitely claim the instability of the vacuum with the second hypercritical magnetic 
field or - which is the same - the instability of such field under the pair creation that 
might provide the mechanism for its diminishing. For the present, we state that the 
first hypervalue 0 is such a value of the magnetic field, the exceeding of which would 
already cause restructuring of the vacuum and demand a profound revision of quantum 
electrodynamics. 

The paper is organized as follows. In Section 2 we revisit Goldstein’s solution by 
referring to various possibilities in approaching the falling to the center, especially the 
one invoked by the previous work [ 13 ], 0 

In Section 3 we derive the ultimate two-dimensional form of the Bethe-Salpeter 
equation in its differential version characteristic of the ladder approximation, when the 
magnetic field tends to infinity, with the help of expansion over the complete set of Ritus 
matrix eigenfunctions m- The latter accumulate the spacial and spinor dependence 
on the transversal-to-the-field degree of freedom. The Fourier-Ritus transform of the 
Bethe-Salpeter amplitude obeys an infinite chain of coupled differential equations that 
decouple in the limit of large B , so that we are left with one closed equation for the 
amplitude component with the Landau quantum numbers of the electron and positron 
both equal to zero, while the components with other values of Landau quantum numbers 
vanish in this limit. The resulting equation is a differential equation with respect to two 
variables that are the differences of the particle coordinates: along the time t = £q — Xq 
and along the magnetic held z = x\ — Xg. It contains only two Dirac matrices 70 and 73 
and can be alternatively written using 2x2 Pauli matrices. Arbitrary external electric 
held E along B is also included, E <C B. By introducing different masses the resulting 
two-dimensional equation is easily modified to cover also the case of an one-electron 
atom in strong magnetic held and/or other pairs of charged particles. 

I 11 Section 4 the ultra-relativistic solutions (possessing maximum symmetry) to 
the equation derived in Section 3 are depicted corresponding to the vanishing energy- 
momentum of the e + e _ -state, and the hrst hypercritical magnetic held is found basing on 
the standing wave boundary condition imposed on the lower border of the normalizing 
volume - as prescribed by the theory in Refs. [ 13 ], 0 Also the standard cut-off 
procedure of Refs. 0 [ 3 ] is fulfilled to give practically the same value 0 . Further, 
we estimate possible modifications that might be introduced by radiative corrections to 
the mass operator, to find that these cannot change the conclusions any essentially, and 
discuss the deconhnement. 
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The fully relativistic Bethe-Salpeter equation for a system of two Fermions of masses 
m a and and opposite charges has in the ladder approximation the form m 

[(i<9a ^a)]a/3[(i^b TXlb X j3v{x a; ^b) ^8>TtOlD rnn ( y X a x b)'~i a g r ) fu/ Xfi 1 ' (. X b-i -^b) • (2) 

Here 7 ’s are the Dirac gamma-matrices, the two-time wave function X( x *> ay,) is a 4 x 4 
matrix with respect to spinor indices, the corresponding Greek letters running the values 
(1,2,3,4). The summation over the repeated vector indices m, n — 0,1, 2, 3 is also meant. 
The derivatives 



7od° b - 


i = 1,2,3 


(3) 


act on the first and the second arguments of X( x a, x b), respectively, a is the fine 
structure constant, and D mn (x a — x b) is the photon propagator. Note that the Feynman 
photon Green function Dp of Ref. iIS] is defined as 2D. The translational invariance 
implies that the solution, which is an eigenfunction of the translation operator, should 
have the form 

Xp( X a, Zb) = exp ^ip Xa + rjp(x), (4) 

where x = x a — Xb- 

Equation © is a complicated set of differential equations, which can, however, be 
essentially simplified, if one assumes = 0 as explained in the previous subsection. In 
this case, equation © can be transcribed as 

(i d - m a )g 0 (x)(-i d - m h ) = -i8TTaD mn (x)'y m ‘r] 0 ( y x )(^ n ) T > (5) 


where the superscript T means transposition and the derivative acts on the relative 
variable x to the right or to the left of it according to what is prescribed by the direction 
of the arrow. With the help of the known relation [3] 

= -C~XnC, ( 6 ) 

where C is the charge conjugation matrix, eq. © is transformed to 

(i d - m a )rjo(x)C~ 1 (i d - m h ) = i87raA nn (x)7 m ?7o(a:)C _1 7 ri . (7) 

The Bethe-Salpeter amplitude rjo{;x)C~ l is the one for the case where the two Fermions 
are Dirac conjugated to one another, i.e. are electron and positron, equation © with 
m a = rrip = m being just the Bethe-Salpeter equation for electron-positron system 4j. 

We take the photon propagator in the Feynman gauge D mn (x) = g mn D(x 2 ), where 
the metric tensor obeys diag g mn = ( 1 ,- 1 ,— 1 ,- 1 ). The Ansatz 

Vo( x ) = ®( x h 5 C, ( 8 ) 

where Q(x) is a unit matrix containing a single function of x, is then consistent with 
the set © and turns it into the equation 

(— □ + m 2 )Q(x) + 32i«7 tD(x 2 )Q(x) = 0, 


(9) 
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where □ = — <9 q + A is the Laplace operator. The bgure 32 here has resulted from the 
multiplication of 8 in © by 4, which is associated with the dimension of the space: 

" = 4 ( 10 ) 


m,n=0,l,2,3 


Qran^mTin 


The photon propagator is singular on the light cone x 2 = 0: 
D(x 2 ) = 


4vr 2 x 2 ' (11) 

In the space-like region x 2 < 0, for the most symmetrical state where the solution does 
not depend on the angles in the 4-dimensional Euclidean space, Q(x) = Q(s), s = v—x 2 , 
equation © with eq. (HU) for D(x 2 ) becomes the Bessel differential equation 

d 2 0 3d© 

s ds 


ds 2 


2 ^ ^ 
m © = —-0. 

7TS 2 


( 12 ) 


Its solution is known as Goldstein’s solution [2] to the Bethe-Salpeter equation with 
P M = 0. (See the review JE], where other ultrarelativistic solutions, corresponding to 
Ansatze different from the above are also listed). Near the singular point s = 0 equation 
ina has the asymptotic form 


d 2 © 3 d© 8cc . . 

ds 2 S ds 7TS 2 " ’ 


(13) 


which is also the asymptotic form of the full (with P M ^ 0) Bethe-Salpeter equation © . 
Its solutions behave near the singular point s = 0 like s a , where 


a= - 1 ± 


8 a 

7T 


(14) 


The substitution 0(s) = T(s)s 2 eliminates the Erst derivative and reduces equation 
m to the standard form, 

d 2 T 1 /8a 3' 

77 


1 /8a 

S* \ 7T 


T = — m 2 4>, 0 < s < 00 , 


(15) 


ds 2 S 2 \ 7T 4 / 

of a Schrodinger-like equation with purely centrifugal - attractive or repulsive, depending 
on the sign of the difference (— — |) - potential. Its solutions behave near the singular 
point s = 0 like 

T(s) ~ (16) 


The same as for the usual radial Schrodinger equation, the natural mathematical 
requirement that the norm f Q |4/(s)| 2 ds be convergent at the lower limit s = 0 is not 
yet sufficient to rule out the more singular solution, which corresponds to the lower 
sign in in the whole range where the square root in Hi, Hi is real (in the 
present case this requirement would separate the less singular solution only for formally 
negative a !). To do so an additional physical requirement concerning the behavior of 
the wave function near the origin is usually imposed to fix the eigenvalue problem in 
quantum mechanics [I[. For the Bethe-Salpeter equation such physical requirement was 
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established by Mandelstam El- It reads that the integral over a small three-dimensional 


closed hypersurface 5(3) of the 4-vector current density 0(x)d /J ,0(x) around the origin 



should tend to zero together with the radius s of the hypersphere 5(3) in the Euclidean 


4-space-time. This implies that 0(s) should increase slower than s 1 and makes only 
solutions with the upper sign acceptable, provided that a < a cr , where 



The above requirement also rules out the both oscillating solutions when a > a cr . If, 
however, one keeps to a weaker condition that the current flow ED be finite as s —» 0 
and, correspondingly, 0 (s) increase no faster than s -1 , the both oscillating solutions are 
compatible with it (note the complex conjugation sign in (fTTI) . due to which the square 
root does not appear in it if a > a CT ). Such situation is typical of the falling to the 
center phenomenon. 

The full Bethe-Salpeter equation © certainly has bound states, when a < a cr , 
corresponding to positronium atom, and the above condition serves to define them. 
The binding energy of the realistic, a = 1/137, positronium makes about a millionth 
fraction of its rest mass. On the contrary, there is no bound state described by equation 
ED- The exact solution to the latter, decreasing at infinity, is expressed in terms of the 
McDonald function as 0(s) = (1 /s)K r —. Its asymptotic behavior near s = 0 

y 7T__ 

is a linear combination of the both terms ED and, hence this solution is forbidden by 
the boundary conditions discussed above, provided that a < a cr . This means that no 
bound state exists with = P = 0 for the coupling, smaller than the critical value 
a = a cr ED, and the gap between electrons and positrons survives down to this value. 

Our main concern is about what happens in the overcritical region of the coupling 
constant a > a cr . Three different approaches have the right to exist for treating this 
case. 

The Erst is the same as the one used for considering a Dirac electron in the Coulomb 
held of a nucleus with its charge greater than 1/a, Z > 137 (see 3], [4]). In that 
approach the finite size of the nucleus was exploited as providing a natural cut-off to 
the singular potential. There is no such natural fundamental length in our problem, but 
if we introduce it formally, for instance by shifting the pole in the photon propagator 
ED from the light cone inwards the time-like domain, D mn (x 2 ) x l/(x 2 — A 2 ) we would 
come to the situation where the positronium gradually approaches the point Po = P = 0 
as a grows and reaches it at certain a cr = a(A). In other words, where there is falling to 
the center, the attraction is so strong that for a sufficient value of the coupling the level 
becomes so deep as to fully compensate for the whole mass of the positronium atom. 
When analogous situation occurs in the above case of the supercharged nucleus, the 
electron level dips into the Dirac sea of positron states, the vacuum becomes unstable 
with respect to e + e _ -pair creation. The essential disadvantage of this approach is that 
all important quantities, the pair production probability among them, depend on the 
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cut-off length and do not have a definite limit when A —> 0. This fact makes the 
results doubtful after the cut-off length becomes less than the characteristic length 
of the problem, which is the electron Compton length m -1 . This is the case for the 
(supercritical) nucleus, whose size is adopted to be 10“ 12 cm -C m~ l . Down to what 
other border should one believe the result once it does not converge? 

The second approach might be dependent on the von Neuman technique of the so 
called self-adjoint extension of the Hamiltonian (see the pioneering works Pi, pi and 
the monograph According to pi and PI there is an infinite number of discrete 

eigenvalues of the operator, whose differential expression is determined by the left-hand 
side of equation (USD- These extend to m 2 = oo for fixed a. The self-adjoint extension is 
defined up to an arbitrary parameter that can always be chosen in such a way as to make 
any given — m 2 - the electron mass squared - an eigenvalue. Then one should conclude 
that the tightly bound states exist beyond the point a = 7 t/8 . The disadvantage of this 
approach is in that no physical criterion is known to fix the arbitrariness of the self- 
adjoint extension p|. If one likes to take the above prescription seriously, one would 
have to alter the choice of this parameter when going to a different a or m 2 . We do 
not know if the method of self-adjoint extension was ever applied to the problem under 
consideration. In the problem of supercritical nucleus the application of this method 
yields the result that the electron level never sinks into the continuum [IH], hence there 
is no place for the electron-positron pairs production effect. 

A third approach is, in our opinion, most adequate. It is to treat equation (USD as an 
eigenvalue problem with respect to the coupling constant a. By putting this eigenvalue 
problem in the form of an integral equation it was demonstrated EH, PI that the 
corresponding integral kernel gives rise to a self-adjoint operator, once an appropriate 
norm is finite, and the eigenvalues a < a cr make a discrete set. If extended to the 
supercritical region a > a cr , this procedure may be thought of to be equivalent to 
the approach recently developed by one of the present authors PI wherein equation 
m is to be represented - by bringing the singular term —8(a/7rs 2 )T to the right-hand 
side, and taking the term —m 2v k to the left-hand side - as the generalized eigenvalue 
problem for the differential operator — ^ + -yA + m 2 defining the spectrum of a. This 
operator is self-adjoint in the (rigged) Hilbert space of functions, orthonorma li zable with 
the singular measure s -2 ds and subjected to the boundary condition 

tf(s 0 ) = 0 (19) 

imposed at the lower edge s = So of a normalization box. As long as Sq is finite we 
face a discrete spectrum of a. The eigen-solutions are standing waves at the lower 
edge of the box and decrease at s —> oo like exp (—ms). In the limit so —*► 0 the 
levels condense to make a continuum of states in the supercritical region a > a CT . The 
norm of the state vector calculated with the singular measure$ diverges in this limit, 

f The requirement of square integrability with the above singular measure, when extended back to the 
values a < a CT (so = 0 in this case), just excludes the lower sign in 1161) and makes the imposing of 
the conditions o discussed above for the bound state problem unnecessary. 
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what makes it possible to normalize the solution to 5-function and hence interpret it as 
corresponding to a free particle living mostly near the singularity. This situation refers 
to the kinematical domain called sector III in m- The corresponding solutions were 
called confined states, since their wave-function decreases at infinity like that of bound 
states. Near the origin there are free particles, incoming from the origin, then scattered 
elastically inwards and then outgoing back to the origin. 

When applied to the supercritical nucleus 5], this approach led to the effect of 
absorption of electrons by the nucleus and to the known effect of electron-positron pair 
production, the corresponding probabilities being calculated in a cut-off-independent 
way. These effects, however, occurred in another kinematical domain, called sector IV, 
where the particles are free also at large distances from the singular center. To see, 
if analogous effects are intrinsic to the positronium in the supercritical case a > a cr 
and cause an instability of the vacuum state relative to the pair production it would 
be necessary to go beyond the kinematical restriction P = 0 , since we need solutions, 
oscillating at infinity as well as near zero to this end. Leaving this task for future, we 
now consider the features of the confined state. 

The solution to equation m for a > a CI that decreases when s —> oo is given by 
the McDonald function with imaginary index 


T(s) = \fs K.^ as—j-(ms). (20) 

This function is real. The lower edge of the normalization box should be chosen much 
smaller than the characteristic length of the problem, which is the electron Compton 
length, So mT 1 . Then one can use the asymptotic form of the McDonald function near 
zero to write the standing wave boundary condition (USD as an equation for determining 
the spectrum of a 


or 


ms 0 

~ 2 ~ 


2v 


r(i + Q 

r*(i - v)’ 



( 21 ) 


v hi —= iarg Y{y + 1) — inn, n — 0, ±1, ±2,... (22) 

Confining ourselves to the values of the coupling constant that do not differ much from 
the critical value, \v\ -C 1 we may exploit the approximation for the Euler T-function 


r(l + v) = 1 — z'C'e, 


(23) 


where Ce = 0.577 is the Euler constant, to get 


In 




n = 1 , 2 ,... 


(24) 


We have expelled the non-positive integers n from here, since they would lead to the 
roots for mso of the order of or much larger than unity in contradiction to the adopted 
condition So m -1 . For such values the asymptotic representation of the McDonald 
function used above is not valid. It may be checked that there are no other zeros of 
McDonald function, besides those in El, enumerated by positive integers. Finally, the 
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discrete spectrum of the coupling constant above the critical value n/8 close to it is 
given as 



(25) 


It is seen explicitly, that the eigenvalues do condense when the lower edge of the 
normalization box s 0 tends to zero. The wave function (ED is mostly concentrated 
inside the Compton radius m -1 , but the probability density is confined to the region 


close to s = s 0 —* 0 due to the singularity of the measure s 2 near this point. 


The solution of the Bethe-Salpeter equation, originally written for positronium 
atom, relates in the ultra-relativistic situation P f/ = 0 considered, as a matter of 
fact to the vacuum as well. Indeed, the state described has no total energy and no 
total momentum and correspondingly does not depend on the coordinate sums of the 
constituting particles. It is not separated from the vacuum by an energy barrier. Besides, 
it is proportional to the unit matrix in the spinor space and is 0(3.1)-invariant, i.e. 
possesses the maximum symmetry, as the vacuum should do. On the other hand, 
this state has a nontrivial space-time structure, described by the dependence on the 
coordinate differences, which implies the concentration of the state near zero value of 
the relative coordinate. These considerations may mean a need of the restructuring of 
the vacuum when the coupling constant exceeds the critical value and serve to establish 
the band of values for this constant beyond which the existing standard concepts no 
longer hold. This is how the things stand with the positronium atom. 

In the next section we derive a two-dimensional analog of equation Q relating to the 
case where a strong magnetic field is imposed, and describe all solutions corresponding 
to the fall-down to the center. The two-dimensioning makes this phenomenon stronger. 
We shall return to the analysis of its consequences in the subsequent section. The 
important difference with the present situation will be that the agent providing the fall- 
down to the center will be the external magnetic field, whereas the coupling constant 
will be kept equal to its experimental value a = 1/137 throughout. 

3. Derivation of two-dimensional Bethe-Salpeter equation in asymptotically 
strong magnetic field 

The view that charged particles in a strong constant magnetic field are confined to the 
lowest Landau level and behave effectively as if they possess only one spacial degree of 
freedom - the one along the magnetic field - is widely accepted. Moreover, a conjecture 
exists [221 that the Feynman rules in the high magnetic field limit may be directly 
served by two-dimensional (one space + one time) form of electron propagators. As 
applied to the Bethe-Salpeter equation, the dimensional reduction in high magnetic 
field was considered in 0,0,1, n. In these references the well-known simultaneous 
approximation to the Bethe-Salpeter equation taken in the integral form was exploited, 
appropriate for nonrclativistic treatment of the relative motion of the two charged 
particles. Once we shall in the next Section be interested in the ultrarelativistic regime, 
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we reject from using this approximation, and find it convenient to deal only with the 
differential form of the Bethe-Salpeter equation. 

The electron-positron bound state is described by the Bethe-Salpeter amplitude 
(wave function) )( a ^(x e , x p ) subject to the fully relativistic equation [TH] . which in the 
ladder approximation in a magnetic field may be written as 

[i<9 e -m + eA(x e )] a0 [id p -nn- eA{x p )\^ X^a; 6 , xP ) = 

= -i8naDij(x e - x p )[^ i ] a p['y j ]„ u x p ) (26) 


Here x e ,x p are the electron and positron 4-coordinates, D i3 (x e — x p ) is the photon 
propagator, and we have explicitly written the spinor indices a, /3, p, u — 1,2, 3,4. 

We refer to, if needed, the so called spinor representation of the Dirac y-matrices 
in the block form 


7o = 


0 

1 



7 k 


0 -cy 

(7 k 0 


(27) 


cy are the Pauli matrices: 


(7 1 = 


0 

1 



kt 2 = 


0 1 

-1 0 



(28) 


k = 1,2,3; 77i is the electron mass, e the absolute value of its charge e = 2-^7ra. 
The metrics in the Minkowsky space is diag g l3 = (1, — 1, — 1, — 1). The vector 
potential of the constant and homogeneous magnetic field £?, directed along the axis 
3 (£>3 = B , £>i i2 = 0), is chosen in the asymmetric gauge 


= -Bx 2 , A), 2,s( x ) = <*. 


(29) 


With this choice, the translational invariance along the directions 0,1,3 holds. 
Solutions to equation (HU) may be represented in the form 


j\j(of,:r p ) = 

V( x o - xg, x‘ s - xg, xg,. xg i2 ) exp{i(P 0 (xg + xg) 


A(xg + xg))}, (30) 


where Po ,3 are the center-of-mass 4-momentum components of the longitudinal motion, 
that expresses the translational invariance along the longitudinal directions (0,3). We 
do not find it convenient to be using explicitly consequences of the magnetic translation 
invariance 22]. (The general representation for the Bethe-Salpeter amplitude that 
incorporates these features may be found in |29, 129). Denoting the differences 
x q — Xq = t, x l — x z = z we obtain the equation 



- m + id e , 

2 - 1 


[rj(t,z, x Y)\p 




A 

e7iHi(a; 2 ) 


—id\\ —-— m + i<9^ + e 7 iHi(x 2 ) 


a/3 

- 


-i8vra D tj (t, z, x e 1>2 - x\ 2 ) [7i] a/ 9 [7jW x Y)\fr ’ 


(31) 


where x± = (x 1 ,x 2 ), -d± = yi<9i + y 2 <9 2 , ((<9j_)i = di, i = 1,2), d\\ = Xyo - <9-73, 
P\\ = PolO ~ ^ 373 - 
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Expand the dependence of solution of equation EB on the transversal degrees of freedom 
into the series over the (complete set of) Ritus [H] matrix eigenfunctions" 1 " Eh(x 2 ) 


W(M,7 P )W = E eipW K-( I 5)];‘K»K)]fe 1 -w [rj h Mt,z)]‘ 


(32) 


h e hp 


Here T)h e h,p(t,z) denote unknown functions that depend on the differences of the 
longitudinal variables, while the Ritus matrix functions e ipiXl Eh(x 2 ) depend on the 
individual coordinates x e { v 2 transversal to the field. The Ritus matrix functions and 
the unknown functions TJh e hp(t , z) are labelled by two pairs h e , h p of quantum numbers 
h = (k,pi, ), each pair relating to one out of the two particles in a magnetic field. The 
Landau quantum number k runs all nonnegative integers k — 0,1, 2, 3..., while pi is the 
particle momentum component along the transversal axis 1. Recall that the potential 
A M (x) flU does not depend on xq, so that pi does conserve. This quantum number is 
connected with the orbit center coordinate x 2 along the axis 2 [T|, pi — —x 2 eB. 

The matrix functions e ipixi Ef p (x 2 ) for transverse motion in the magnetic field (I2H1) . 
relating in © to electrons (e) and positrons (p), are 4x4 matrices, formed, in the 
spinor representation, by four eigen-bispinors of the operator (i<9j_ ± eA) 2 

(A ± eA)%,e' p ™ R p (x 2 )]7> = —2eBke ipixl {E‘»(x 2 )]<ry (33) 

placed, as columns, side by side [13J. Here the upper and lower signs relate to electron 
and positron, respectively, while a = ±1 and 7 = ±1 are eigenvalues of the operators 


S 3 


03 0 \ 

0 (x 3 y’ 



(34) 


diagonal in the spinor representation, to which the same 4-spinors are eigen-bispinors'* 


i-rf 7 ’ = 7 E { r\ 




3 




= crE, 


0,7) 


sg 


(35) 


The couple of indices a = (<7,7) is united into one index a in the expansion (|32|) . 
a = 1,2, 3,4 according to the convention: (+1,-1) = 1, (-1,-1) = 2, (+1,+1) = 

3, (—1,+1) = 4. With this convention, the set of 4-spinors [E , /l (x 2 ]| t <T ’ l) = E^(x 2 )“ can 
be dealt with as a 4 x 4 matrix, the united index a spanning a matrix space, where the 
usual algebra of 7 -matrices may act. Correspondingly, in (fT^l) the unknown function 
[kjh B hp{t, £)]“ e “ P is a matrix in the same space, and contracts with the Ritus matrix 
function. 

Following [S|, the matrix functions in expansion (H?!) can be written in the block 
form as diagonal matrices 


gipi^l 


E e f P (x 2 ) 


a e ’P(/r;x li2 ) 0 

0 a e,p (h-x i j2 ) 




+ Our definition of the matrix eigenfunctions differs from that of Ref. d in that the longitudinal 
degrees of freedom are not included and the factor exp(ipia:i) is separated. 

* Henceforth, if superscripts e or p are omitted, the corresponding equations relate both to electrons 
and positrons in a form-invariant way. 
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i e iP 


(h',x i )2 ) = 


(36) 


a+i(/i;x i, 2 ) 0 

0 a!^(/i;xi )2 ) 

Here a® ,p (/r; a; li2 ) are eigenfunctions of the two (for each sign ±) operators 
[((i<9_i_)a ± eA a ) 2 =f creB\, labelled by the two values a = 1, — 1 


[- ((id±)a ± eA Q ) 2 ± <rei?]<£ p (/i; xi, 2 ) = -2eBka e j p (h; xi :2 ), 


(37) 


Namely, (we omit the subscript ”1” by p\ in what follows) 

«r(h; x 1)2 ) = e ipxi U k+ ±^ [VAb (x 2 ±A)), k = 0,1, 2..., (38) 

with 

U n (0 = exp |-|| (2 n n!V5F)-^ n (0 (39) 

being the normalized Hermite functions (H n (£) are the Hermite polynomials). Equations 
(EH) are the same as El due to the relation 

(id ± ± eAf = - ((i<9j_) a ± eA a f ± eHE 3 (40) 


and to eq. dSSD- Besides, the matrix functions m are eigenfunctions to the operator 
—\d\ , as commuting with S 3 and 75 El, and with (i<9j_ + eA ) 2 . The corresponding 
eigenvalue pi does not, however, appear in the r.-h. side of El due to the well-known 
degeneracy of electron spectrum in a constant magnetic field. 

The orthonormality relation for the Hermite functions 

/ OO 

U n {t)U n '(Z)d£ = 8 nn '. (41) 

-00 

implies the orthogonality of the Ritus matrix eigenfunctions in the form 

V^B J E* h (x 2 )*E h '(x 2 )idx2 = 5 kk ,5 aa (42) 

As a matter of fact, the matrix functions Eh(x 2 ) are real, and we henceforth omit the 
complex conjugation sign ”*”. 

The matrix functions e ipxi E^ p (x 2 ) El commute with the longitudinal part 
± i9|| — Py/2 — m of the Dirac operator in (131 H . owing to the commutativity property 

[E h (x 2 ), 7o, 3 ]-=0, (43) 

and are ra. in a sense, matrix eigenfunctions of the transversal part of Dirac operator 
(not only of its square El) 

(i<9j_ ± eA)e ipxi E e *(x 2 ) = ±V2eBk e ipxi E e k P (x 2 )-f 1 . (44) 

The Landau quantum number k appears here as a ’’universal eigenvalue” thanks to 
the mechanism, easy to trace, according to which the differential operator in the left- 
hand side of eq. El acts as a lowering or rising operator on the functions m, whereas 
the matrix cr 2 , involved in q 2 , interchanges the places the functions U k , £4_i occupy 
in the columns. Contrary to relations, which explicitly include the variable cr, whose 
value forms the number of the corresponding column, relations El, El, El, and 
the first relation in El are covariant with respect to passing to other representation of 
7-matrices, where the matrix Eh(x 2 ) may become non-diagonal. 
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Now we are in a position to use expansion (HH in equation dm We left multiply it 
by (27 t )~ 2 eB e~ lpax i Ef- e (x 2 )e~ ipPxP E2 P (x 2 ), then integrate over d 2 x ® 2 d?x\ 2 . After using 
(1-H1) and (14511 . and exploiting the orthonormality relation (14211 for the summation over 
the quantum numbers h e,p = (k e,p , pV), the following expression : 




m + 71 V 2 eBk e 


J aot e 


—i9|| —— m — 71 V2 eBkP 


[Tj^hp(t,z)} a 01 ( 45 ) 


HaP 


is obtained for the left-hand side of the Fourier-Ritus-transformed equation (ETTll . We 
omitted the bars over the quantum numbers. 

Taking the expression 


Dij [t,z,x\ 2 -x\ 2 ) = t|^( t 2 


Xi 


X 


Pl 2 


Xn 


X 


P\2 


( 46 ) 


a 

2n 3 


for the photon propagator in the Feynman gauge, we may then write the right-hand 
side of Ritus-transformed equation m as 

[ dp e dp p 9ij f [E^{x e 2 )zfiE^(x e 2 )\ aae [^(a^bi-^O^W [T] hehP (t, z)] aBaP 

^ k e kP ^ 

e i{ P e -F)xi e i(p p -F)^i e Pd 2 x® 2 d 2 x p l 2 


Z 2 + (x® — x\) 2 + (x| — X 2) 2 — t 2 


( 47 ) 


Integrating explicitly the exponentials in (El) over the variable A" = (x® + x^)/ 2 , we 
obtain the following expression: 


7 r 


f dp dP 1 HP ! - PHJ2 3<s f I E^(x‘ 2 HE‘ h .(4)U. [£&(*?) 'fjEU4)U- ■ 

•* k e kP J 

exp(ix(p — p))dx 


[ 7 /h e hp(T Z )Y 


eRdxgdxg, 


( 48 ) 


Z 2 + X 2 + (x| — X 2) 2 — t 2 

where the new integration variables x = x® — x p , Pi = p e + p p , p = (p® — p p )/2 and 
the new definitions Pi = p® + p p , p = (p e — p p )/2 have been introduced. The pairs of 
quantum numbers in El) are 


r’ p = (r’ p ,P±p), 

Hence the arguments of the functions 


A*' P =(V’ p .y±p). 


Ill 


are: 


VeB 


Xo T 


Pi + 2 p 

2eB 


VeB \ 




x% + 


Pi + 2 p 

2eB 




VeB 


Xn 


Pi - 2 p 

2eB 


, VeB (x 2 




( 49 ) 


Pi ~ 2 p 

2eB 


successively as the functions P^(x i i2 ) appear in (14811 from left to right. After fulfilling 
the integration over dPi with the use of the 5-function, introduce the new integration 
variable q = p — p instead of p, and the integration variables xfj = x 2 + (Pi + 2p)/2eB, 
x 2 = x 2 — (Pi — 2p)/2eB instead of x| and x p . Then ea. dlHll may be written as 
a 


n 


/d q 9ij [T)h*hp(t,z)Y 

^ k e kP ^ 


( 50 ) 
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exp(—i xq) dx eB dx | dx p 


z 2 + x 2 + [x‘ 


—p _ Pi-q 

eB 


x 


-1 2 


Now the pairs of quantum numbers in (Eu are 


5*' p = (F' P , P±v), 


Pi 


IP P = (fc e ’ p , -y ±g±P). 


Hence the arguments of the functions (|HHD in (ED from left to right 


are: 


VeBx 2, Ai/eH , \/eH ("xl)-, 


eH 


eH 


VeBx 2 


(51) 


(52) 

(53) 


5.5. Adiabatic approximation. 

Now we aim at passing to the large magnetic field regime in the Bethe-Salpeter equation, 
with (ED as the left-hand side and ED as the right-hand side. Define the dimensionless 
integration variables w = x\feB , q 1 = q/VeB, £ e,p = x% p y/eB in eq. ED- Then it takes 
the form 

£ /<* £ * 

J k e kP 



exTp(—iwq')dwd^ e d^ p 


z 2 + 


eB ' eB ^ 


e p 


Pi 

\feB 


9 



The pairs of quantum numbers in ED are 

v .p = (F ,p j Hip), ft ..p = (*«>, b- ± c/VTb ±p). 
The arguments of the functions ED in ED from left to right are: 

f, f + V e-q', e- 


(54) 


(55) 

(56) 


When considering the large field behavior we admit for completeness that the 
difference between the centers of orbits along the axis 2 X 2 — x p = —^ may be 
kept finite, in other words that the transversal momentum P\ grows linearly with the 
field. We shall see that that big transversal momenta do not contradict dimensional 
compactikcation, but produce an extra regularization of the light-cone singularity. 

In the region, where the 2-interval (z 2 — f 2 ) 1 / 2 essentially exceeds the Larmour radius 
Lb = 1/VeB, 


z 2 - t 2 » L\ (57) 

one may neglect the dependence on the integration variables w and later on £ 6jP in the 
denominator. Integration over w produces 27 t 5(</), which annihilates the dependence 
on q' in the arguments ED of the Hermite functions, and they all equalize. 

Let us depict this mechanism in more detail. Fulfill explicitly the integration over 
dw in ED : 


exp(— iwq') die 
.2 _ +2 1 ni 4. df 

b ' eB ' eB 


V eB 


IT 


~t 2 + 


eB 


(o(q) exp (- q'VeB(z 2 - t 2 ) + A 2 ) + 
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O(-q') exp ( q f y/eB(z 2 - t 2 ) + A 2 ) j , (58) 

where 

a 2 = fe -e- 

and d(q') is the step function, 

{ 1 when q' > 0, 

| when q 1 = 0, (60) 

0 when q' < 0. 

Due to the decreasing exponential in (El the variables £ e,p do not exceed unity in the 
order of magnitude and can be neglected as compared to ^== in El- Unless q 1 is large 

p 2 

it may be neglected as compared to the same term in El, too. Then A 2 = and after 
El is substituted in El and integrated over d q' the contribution comes only from the 
integration within the shrinking region \q'\ < ( eB[z 2 — t 2 + ])~ 1 • Then q' can be also 

neglected in the arguments El- If, contrary to the previous assumption, we admit that 
\q'\ is of the order of -j= ~ \JeB we see that the exponentials in El fast decrease with 
the growth of the magnetic field as exp(— eB(z 2 — f 2 )), and therefore such values of \q'\ 
do not contribute to the integration. If we admit, last, that \q'\ |4=|, we find that 

the contribution exp(— \q'\y/eB(z 2 — t 2 ) + (q') 2 ) from the integration over such values 
is still smaller. Thus, we have justified the possibility to omit the dependence on q' in 
El and in El, and also on £ e,p in El- Now we can perform the integration over d q' 
to obtain the following expression for El 
2«7r _1 


Pi 

y/eB 


(59) 


+ 


p? 

(eP)2 


t 2 


Y,g,. m-(^) 7^.te)U<r f6i) 


fc e fcp 


It remains yet to argue that the limit El is valid also when the term j s no t kept. 
In this case we no longer can disregard q' inside A 2 when q' is less than or of the order 
of unity. But we can disregard A 2 as compared with eB(z 2 — t 2 ) to make sure that the 
integration over d q' is restricted to the region close to zero \q\' < ( eB(z 2 — t 2 )) -1 / 2 and 
hence set q' = 0 in El- The contribution of large q' is small as before. 

The integration over £ 6jP of the terms with i = 0,3 in El yields the Kroneker 
deltas S ke j ^ 5 k pjp due to the orthonormality El of the Hermite functions thanks to the 
commutativity El of the Ritus matrix functions El with 7 q and 73. On the contrary, 
71, 72 do not commute with El- This implies the appearance of terms, non-diagonal in 
Landau quantum numbers, like 5 ke -^ ±1 and S kPk p ±1 , in El , proportional to (i — 1,2) : 

r£ ±1 ,p. ±I = £ f [^(4bi^.(4)Udf 

fc e fcP ^ 


[Epp( x 2 ) 1 iE kP (^2 )]d£ p [Tjh*hp(t, z)] 


aa y _ 


= £ 

k e kP 


0 - A r,e \ 

A Ke 0 ) 


Ai P 


—A* 

0 





aa e 


(62) 











Positronium collapse 


18 


Here x 2 P are expressed in terms of £ through the chain of the changes of variables 
made above starting from ED, so that all the arguments of the Hermite functions have 
become equal to £. Besides, 

h e ’ p = (fc e ’ p ,p e ’ P ), 7T’ P = (F’ P ,F’ P ), P e + P p = Pi. (63) 


= / a'(h,x 2 ) a'{h,x 2 )d£, i = 1,2 


(64) 


Ai 1} = 

kk 


0 a' +1 (h, x 2 )a / _ 1 (/c,a:2) 

a'_ x {h,X 2 )a' +1 (h,x 2 ) 0 

0 

% fc+i ^ 


d£ = 


"fc, fc-i 


(65) 




0 —a( fl (/i, x 2 )a / _ 1 (fc, x 2 ) 

X 2 )a' +1 (h,x 2 ) 0 

0 -St 


d£ = 


= i 


5: 


k, k-\-l 


J k, k—1 

o 


( 66 ) 


The prime over a indicates that the exponential exp(ipxi) is dropped from the definitions 
(ED and HD- The non-diagonal Kronecker deltas appeared, because a' ±1 (h,x 2 ) are 
multiplied by a' Tl (h,x 2 ) under the action of the er^-blocks in 7 1j2 (ED- In the final 
form, the matrices in (E2D are 


0 

Ai 

kk 


Q U ) — 9 (7l(d=<5fc ) fc_i + $k,k+l) + ^k,k+l)) , 


(67) 


with the upper sign relating to i = 1 and the lower one to i = 2. Now equation m 
acquires the following form, 


id|l —— m -f 71 V 2 eBk e 


J OiOT L 


—i<9|| —di — m — 71 F 2 eBk p 


[Plh e hp(ty ^)] 


_ 


J /ia p 


2«7T 


-1 


p2 

z 2 H- L 


t 2 


E s«n aa e [di]/iQ: p [7//i e /iP (tj ^)] 


^ fcp ±1 ) , P1+P P ! = 

1 = 1,2 


(eB) 2 ^ V=0,3 

The bars over quantum numbers are omitted. This equation is degenerate with respect 
to the difference of the electron and positron momentum components p = (p e — p p )/2 
across the magnetic field, but does depend on its transversal center-of-mass momentum 
Pi = (jf-\-p p ). This dependence is present, however, only for sufficiently large transverse 
momenta P\. 

At the present step of adiabatic approximation we have come, for high magnetic 
field, to the chain of equations (168(1 . in which the unknown function for a given pair 
of Landau quantum numbers k e , k p is tangled with the same function with the Landau 
quantum numbers both shifted by ±1 (in contrast to the general case of a moderate 
magnetic field, where these numbers may be shifted by all positive and negative integers). 
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To be more precise, the chain consists of two mutually disentangled sub-chains. The 
first one includes all functions with the Landau quantum numbers k e , k p both even 
or both odd, and the second includes their even-odd and odd-even combinations. We 
discuss the first sub-chain since it contains the lowest function with k e = k p = 0. We 
argue now that there exists a solution to the first sub-chain of equations (ISED, for which 
all Tj k e, p e. k p p p(t, z) disappear if at least one of the quantum numbers k e ,k p is different 
from zero. Indeed, for k e = k p = 0 equation (jfiHl) then reduces to the closed set 

[i<9|| - y - rn\ aa e [-idy - y - m] MaP [Pjo^o^t, z)]" e " P = 

- 9iih] a ^h]^ [ r ]o,pt;o,p p 1 (t,z)] ae °‘ P , p\+p\ = Pi- (69) 

+ *=0,3 

In writing it we have returned to the initial designation of the electron and positron 
transverse momenta p e { v . Denote for simplicity TJ k e k p = ^7fc e ,pf;fc P ,pf (t, z )- If we consider 
equation (16811 with k e = k p = 1 for Tj n we shall have a nonzero contribution in the right- 
hand side, proportional to TJoo coming from T k . e _ lkP _ 1 , since the other contributions 
k]n,Tf 22 ,P/ 20 , P]o 2 are vanishing according to the assumption. As the left-hand side of 
equation (jSHD now contains a term, infinitely growing with the magnetic held B , it can 
be only satisfied with the function Tjn, infinitely diminishing with B in the domain (11771) 
as 

[P]n] ai “ =~2^b 2 a p/ -- E^ii[7i7i]aa-[7i7i]^ [Tjoo]^ , (70) 

Z + jpWp t *=0,3 

in accord with the assumption made. Thus, the assumption that all Bethe-Salpeter 
amplitudes with nonzero Landau quantum numbers are zero in the large-held case is 
consistent. We state that a solution to the closed set (EH) for 77o, p po,pP (t, z) with all the 
other components equal to zero is a solution to the whole chain (EH). 

The derivation given in this Subsection realizes formally the known heuristic 
argument that, for high magnetic held, the spacing between Landau levels is very large 
and hence the particles taken in the lowest Landau state remain in it. Effectively, 
only the longitudinal degree of freedom survives for large B , the space-time reduction 
taking place. Equation (EH) is a fully relativistic two-dimensional set of equations 
with two space-time arguments t and z and two gamma-matrices 70 and 73 involved. 
Since, unlike the previous works 0 , q, nu, neither the famous equal-time Ansatz 
for the Bethe-Salpeter amplitude ra, nor any other assumption concerning the 11011- 
relativistic character of the internal motion inside the positronium atom was made, the 
equation derived is valid for arbitrary strong binding. It will be analyzed for the extreme 
relativistic case in the next Subsection. 

The two-dimensional equation (1(711) is valid in the space-like domain (ElD- it 
is meaningful provided that its solution is concentrated in this domain. In 11011- 
relativistic or semi-relativistic consideration it is often accepted that the wave function is 
concentrated within the Bohr radius a 0 = (am)^ 1 = 0.5 x 10 -8 cm. It is then estimated 
that the corresponding analog of asymptotic equation (R71 holds true when ao Lb, 
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i.e. for the magnetic fields larger than a 2 m 2 /e ~ 2.35 x 10 9 Gauss. This estimate, 
however, cannot be universal and may be applicable at the most to the magnetic fields 
close to the lower bound where the value of the Bohr radius can be borrowed from the 
theory without the magnetic held. Generally, the question, where the wave function is 
concentrated, should be answered a posteriori by inspecting a solution to equation m- 
Therefore, one can state, how large the fields should be in order that the asymptotic 
equation ( 1 HHI) might be trusted, no sooner that its solution is investigated. We shall 
return to this point when we deal with the ultra-relativistic situation. 

Remind that the transverse total momentum component of the positronium system 
is connected with the separation between the centers of orbits of the electron and 
positron P 1 /(eB) = x\ — x 1 ^ i n the transversal plane, so that the ’’potential” factor 
in eq. ( 1691 ) may be expressed in the following interesting form 

_“_ ( 71 ) 

(r,.e _ T PU _ (, T e _ ,-v.P^2 _ (™e _ ~ P\2’ ' ' 

\ x 0 - L 0J \ x 3 x 3/ V x 2 x 2) 

(cf the corresponding form of the Coulomb potential in the semi-relativistic treatment of 
the Bethe-Salpeter equation in 0, P, m- the difference between the potentials in 0, 
0 and PH lies within the accuracy of the adiabatic approximation). The appearance 
of P f in the potential determines the energy spectrum dependence upon the momentum 
of motion of the two-particle system across the magnetic held like in 0, i, m, M- 
We shall need equation m in a more convenient form. First, transcribe it as 


^ P\\ tr- Pu 

( i<9 ll “ Y - m) Tj 0 j) e . 0 p p(t, z) (—i 9 || - — - m) T = 

=- Vojiiotffaz) t p e i+P V i=Pi- ( 72 ) 

Z +^W~ t '=0,3 

Here the superscript T denotes the transposition. With the help of the relation 
7 J = —C~ lr jiC, with C being the charge conjugation matrix, C 2 = 1, and the anti¬ 
commutation relation [71,75]+ = 0 , 75 = — 1 , we may write for a new Bethe-Salpeter 
amplitude 0 (t,z), defined as 


Q(t,z)= rjo,pt; 0 ,pl(t,z) G75, ( 73 ) 

the equation 

(i^ll - y - m) 0 (t, ^)(—i 5 || - y -m) = 

= -^2- 9 anfi®(t, zhi, P e i+Pi = Pi■ ( 74 ) 

- 2 + w -^ 2 

The unknown function 0 here is a 4 x 4 matrix, which contains as a matter of fact 
only four independent components. In order to correspondingly reduce the number of 
equations in the set (ED, one should note that the 7-matrix algebra in two-dimensional 
space-time should have only four basic elements. In accordance with this fact, only 
the matrices 70,3 are involved in ED- Together with the matrix 7073 and the unit 
matrix I they form the basis, since 70,3 ■ 7o73 = 73,0, 7o = “73 = (7o73) 2 = 1, 
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[70,73] + = [ 7 o,3) 7073]+ = 0 . Using this algebra and the general representation for 
the solution 

0 = al + 670 + C73 + d7o73, ( 75 ) 

one readily obtains a closed set of four first-order differential equations for the four 
functions a, b, c, d of t and z. The same set will be obtained, if one replaces in eq. ED 
and ED the 4 x 4 matrices by the Pauli matrices ED, subject to the same algebraic 
relations, according, for instance, to the rule: 7o =>■ 03, 73 => icr 2 , 7073 =>■ ay. Then 
equation ED becomes a matrix equation 

(idtV 3 + dtcr 2 - y<t 3 + yicr 2 - m)d(t, z)(-iSo- 3 - c^a 2 - y<7 3 + yia 2 - m) = 
2q(7t — ^ 

=-{^30(*, 2)^3 + ^)cr 2 }, Pi + Pi = Pi. (76) 

- 2 +w 

for a 2x2 matrix 1? 

1? = a/ + ba 3 + ica 2 + day. ( 77 ) 

Here / is the 2 x 2 unit matrix, and functions a, b, c, d are the same as in ED- 

The following remark is in order. The derivation above relates to the system of 
two spinor fields, marked by the superscripts e and p, with opposite charges ±e and, 
generally, different masses. Although we kept the masses equal above, it is easy to 
restore their difference by setting m = m e in the left differential operator and m = m p 
in the right one starting with eq. ED throughout. The corresponding Bethe-Salpeter 
amplitude T] is the translationally invariant part ( 1301 ) of the matrix element of the 
chronological product of the spinor field operators 

*>(*•, * p ) =<0|T(^(V)«(i»))|P> ( 78 ) 

between the vacuum < 0 | and the bound state |P >. Once we restrict ourselves 
to the case where one of the particles, "0 e , is an electron and the other, -0 P , is a 
positron, we should take V ;P — Cil> in ED and keep the masses equal. Then the 
Bethe-Salpeter amplitude of two arbitrary opposite-charged fermions y and the electron- 
positron Bethe-Salpeter amplitude p =< O|T(' 0 ^(:r e )' 0 ®(a; p ))|P > are connected as 
q C T = — g C. It follows from (T7H1) that the translationally-invariant part of the 
Ritus transform of g obeys the same equation as ED, but with the sign in front of the 
hatted terms in the right-most Dirac operator reversed, as well as the common sign in 
the right-hand side. The subsequent 75-transformation in ED is useful, since it gives 
the possibility to form the Laplace operator in the subsequent equations. 

3-4- Including an external electric field 

Let us generalize the two-dimensional Bethe-Salpeter equation obtained in the presence 
of a strong magnetic field by including an external electric held, parallel to it, that is 
not supposed to be strong, E -C B. To this end we supplement the potential ED in 
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equation m by two more nonzero components 
A a (x 0 , x 3 ), A 3 (i 0 , Xz)) # 0, 


(79) 


that carry the electric field - not necessarily constant - directed along the axis 3. We 
shall use the collective notations A\\ = (A 0 ,A 3 ), x\\ = {xq,x 3 ), <9jj’ p = <9 q P 7 o — ^ 3 P 73, 
A|| = A 0 7 o — ^373- We shall not exploit now a representation like (EHD, but deal 
directly with the Bethe-Salpeter amplitude )((x e ,x p ) as a function of the electron and 
positron coordinates, and with its Fourier-Ritus transform )(h e hp( x jj,^) connected with 
ij 1 ; Xj_, x v L ) in the same way as (1321) . In place of equation (HTTll one should write 


i9|| C — eA\\{xf) — 7 n + id c ± — e^Ai{x%) i<9|| P + eA^(x?,) — m + i<9^ + e^A^x^) 

J a(3 L 


fllS 


X( x l ’V± P )W = -&™D i:j (t,z,xl 2 -xl 2 ) [7i]o^ [7 jV [X(zJ’ p ,^ p )W , (8°) 

Thanks to the commutativity (TEA the rest of the procedure of the previous Subsection 
remains essentially the same, and we come, in place of m, to the following two- 
dimensional equation 

[Xo,pf;0,Pj (R’iI 5 % || )] 


/s 0 /\ 

idn — eA\\(xf\) — m 

[i<9i| P + eA\\ fx?) — m 

. 

afi L 11 11 11 J 




2q.it 


-i 


z 2 + 


t 2 £ 


9a [7 i\a/3 [Tj'W [Xo,pf ;0,pj ( X *\\ ’ ^f) 


$v 1 


(81) 


(eSR _ L *=0,3 

for a positronium atom in a strong magnetic field placed in a moderate electric field, 
parallel to the magnetic one. In order to apply this equation to a system of two different 
oppositely charged particles interacting with each other through the photon exchange 
and placed into the combination of a strong magnetic and an electric field in the same 
direction, say a relativistic hydrogen atom, one should only distinguish the two masses 
in the first and second square brackets in the left-hand side. 


4. Ultra-relativistic regime in a magnetic field 


In the ultra-relativistic limit, where the positronium mass is completely compensated by 
the mass defect, Pq = 0, for the positronium at rest along the direction of the magnetic 
field P 3 = 0, the most general relativistic-covariant form of the solution m is 

0 = /$ + 9||4 ) 2 + 7o73$3- (82) 


The point is that 7073 is invariant under the Lorentz rotations in the plane (t,z). 
Substituting this into (d with P 0 = P 3 = 0 we get a separate equation for the singlet 
component of (1H21) 


(-□2 + m 2 ) $(t,z) 


4an 1 <!>(t,z) 



-t 2 


and the set of equations 


(□2 + 77i 2 ) $ 3 (t,z) 


4ct7T 1 <4> 3 (t, ^) 


z 2 + 


p? 

(Tup 


t 2 


J 


(83) 
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(-□2 + m 2 )d t § 2 + 2mi<9 z $ 3 = 0, 

(-□2 + m 2 )d z <& 2 + 2mid t $ 3 = 0 (84) 


for the other two components. Here IZb = —d 2 /dt 2 + d 2 jdz 2 is the Laplace operator in 
two dimensions. Note the ’’tachyonic” sign in front of it in the first equation 

Let us differentiate the second equation in (IHUl over z and the third one over t and 
subtract the results from each other. In this way we get that □2^*3 = 0. This, however, 
contradicts the first equation in PD- Therefore, only $3 = 0 is possible. Then, the 
two second equations in (1H31) are satisfied, provided that (— D 2 + m 2 ) < h 2 = 0. We shall 
concentrate in equation Ann in what follows. 

The longitudinal momentum along aq, or the distance between the orbit centers 
along X 2 , plays the role of the effective photon mass and a singular potential regularizator 
in equation PD- The lowest state corresponds to the zero value of the transverse total 
momentum = 0. In this case the equation for the Ritus transform of the Bethe- 
Salpeter amplitude finally becomes 


(-□2 + m 2 ) <L(f, z) 


4 z ) 

7 t(z 2 — t 2 ) 


(85) 


4-1. Fall-down onto the center in the Bethe-Salpeter amplitude for high magnetic field. 
First hypercritical field. 

We are going now to consider the consequences of the fall-down onto the center 
phenomenon present in equation (USD, formally valid for an infinite magnetic field, and 
the alterations introduced by its finiteness. 

I11 the most symmetrical case, when the wave function <h(a;) = < L(s) does not depend 
on the hyperbolic angle 0 in the space-like region of the two-dimensional Minkowsky 
space, t = s sinh 0, z = s cosh 0, s = v z 2 — t 2 equation (|HSD becomes the Bessel 
differential equation 

d 2< 3> 1 d<h 2 4 a 

- — -— + m 2 <f> = —-<L. (86) 

ds^ s ds tis z 

It follows from the derivation procedure in the previous Section 3 that this equation is 
valid within the interval 
1 


VeB 


•C so < s < 00, 


(87) 


where the lower bound Sq depends on the external magnetic held - it should be larger 
than the Larmour radius Lb = ( eB )~ 1 / 2 and tend to zero together with it, as the 
magnetic held tends to infinity. The stronger the held, the ampler the interval of validity, 
the closer to the origin s = 0 the interval of validity of this equation extends. If the 
magnetic held is not sufficiently strong, the lower bound So falls beyond the region where 
the solution is mostly concentrated and the limiting form of the Bethe-Salpeter equation 
becomes noneffective, since it only relates to the asymptotic (large s) region, while the 
rest of the s-axis is served by more complicated initial Bethe-Salpeter equation, not 
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reducible to the two-dimensional form there. This is how the strength of the magnetic 
held participates - note, that the coefficients of ea. (j%Hl) do not contain it. 

In treating the falling to the center below we shall be using sq as the lower edge of 
the normalization box (see the discussion in Section 2). For doing this it is necessary that 
So be much smaller than the electron Compton length, the only dimensional parameter 
in equation m- In this case the asymptotic regime of small distances is achieved 
and nothing in the region s < So beyond the normalization volume - where the two- 
dimensional equations m, d, d, d and hence d are not valid- may affect the 
problem, because this is left behind the event horizon. 

In alternative to this, we might treat So as the cut-off parameter. In this case we 
have had to extend equation m continuously to the region 0 < s < s 0 , simultaneously 
replacing the singularity s~ 2 in it by a model function of s, nonsingular in the origin, 
say, a constant Sq 2 . In this approach the results are dependent on the choice of the 
model function which is intended to substitute for the lack of a treatable equation in 
that region. Besides, the limit so —*► 0 does not exist. The latter fact implies that the 
approach should become invalid for sufficiently small s 0 , i.e., large B. We, nevertheless, 
shall also test the consequences of this approach later in this section to make sure that 
in our special problem the result is not affected any essentially. 

The crucial difference of (1%HI) from equation m, which relates to the case where the 
magnetic held is absent, is the coefficient 1 in front of the first-derivative term instead 
of 3 (this coefficient is D — 1, where D is the dimension of the Minkowsky space in the 
radial part of the Laplacian Do = s~ D+1 d/ds(s D ~ 1 d/ds)). The coefficient 4 in front of 
a in d and (Hi instead of 8 in m is also of geometric origin: the relation 



( 88 ) 


was used when we passed from to (IHH1) , whereas relation m was exploited to pass 
from (J2J to (I12jh Solutions of flHHD behave near the singular point s = 0 like s a , where 



(89) 


The fall-down onto the center [.fj occurs, if a > a cr = 0, i.e., unlike (TOfll . for 
arbitrary small attraction, the genuine value a = 1/137 included. With the substitution 
4>(s) = ^(s)/y/s equation (jHED acquires the standard form of a Schrodinger equation 



so < s < oo, s 0 » (eB) 1 ^ 2 (90) 


Treating the applicability boundary So of this equation as the lower edge of the 
normalization box, as discussed above, s 0 m _1 , we impose the standing wave 
boundary condition d on the solution of d 



( 91 ) 


that decreases at infinity. 
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Starting with a certain small value of the argument ms, the McDonald function 
with imaginary index m oscillates, as s —> 0, passing the zero value infinitely many 
times. Therefore, if s o is sufficiently small the standing wave boundary condition m 
can be definitely satisfied. Keeping to the genuine value of the coupling constant 
a = 1/137 (y = 0.096 i) one may ask: what is the largest possible value s™ ax of s 0 , 
for which the boundary problem ©, m can be solved? By demanding, in accord 
with the validity condition (EZ1) of equation EHD, EHD, that the value of s“ ax should 
exceed the Larmour radius 

(e-B)' 1/2 (92) 


> 


e O^) 2 


one establishes, how large the magnetic held should be in order that the boundary 
problem might have a solution, in other words, that the point Pq = P = 0 might belong 
to the spectrum of bound states of the Bethe-Salpeter equation in its initial form (HID. 

One can use the asymptotic form of the McDonald function near zero to see that 
the boundary condition m is satisfied provided that 


pms o 
V 2 


2u 


r(i + Q 

r*(i — v) 


or 


ms 0 


(93) 

(94) 

Once \v\ is small we may exploit the approximation (fTT!) for the T-function to get 

n — 1,2,... (95) 


z/ln-= iarg Y{y + 1) — inn, n = 0, ±1, ±2, 




a 


We have expelled the non-positive integers n from here for the same reasons as in Section 
2 (see eq.E3D). The maximum value for so is provided by n = 1. The Euler constant 
Ce=0.577 contribution is small as compared to (l/2)yV 3 /a = 32.588. We finally get 


In 


msn 





2 V "a 


or 


So liax = — exp 
m 



n u 

a 


a 


~ — exp{—33} ~ 10 14 —. 
m m 


(96) 


This is fourteen orders of magnitude smaller than the Compton length m 1 = 3.9 • 
10 _11 cm and makes about 10~ 25 cm. Now, in accord with E2D, if the magnetic held 
exceeds the first hypercritical value of 

5 hpcr = ex P + 2C E j ^ 1.6 x 10 28 Bo, (97) 

the positronium ground state with the center-of-mass 4-momentum equal to zero 
appears. Here Bq = m 2 /e = 2.17 x 10 21 cm~ 2 is the Schwinger critical held, or 
B 0 = m 2 c 3 /eh = 1.22 x 10 13 Heaviside-Lorentz units, or B 0 = 4.4 x 10 13 Gauss 
(a = e 2 /Anhc, e = 4.8-10 _1 o a/47t CGSE). Excited positronium states may also reach the 
spectral point P = 0, but this occurs for magnetic helds, tens orders of magnitude larger 
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than © - to be found in the same way from (ESI) with n = 2, 3... The ultra-relativistic 
state P /t = 0 has the internal structure of a confined state, i.e. the one whose wave 
function behaves as a standing wave combination of free particles near the lower edge of 
the normalization box and decreases as exp (—ms) at large distances. The effective ’’Bohr 
radius”, i.e. the value of s that provides the maximum to the wave function m makes 
«max = 0.17m -1 (this fact is established by numerical analysis). This is certainly much 
less than the standard Bohr radius (e 2 m) -1 . Taken at the level of 1/2 of its maximum 
value, the wave-function is concentrated within the limits 0.006 m -1 < s < 1.1 m~ l . 
But the effective region occupied by the confined state is still much closer to s = 0. 
The point is that the probability density of the confined state is the wave function 
squared weighted with the measure s -2 ds singular in the origin m, n and is hence 
concentrated near the edge of the normalization box So = 10 -25 cm, and not in the 
vicinity of the maximum of the wave function. The electric fields at such distances are 
about 10 43 volt/cm. Certainly, there is no evidence that the standard quantum theory 
should be valid under such conditions. This remark gives the freedom of applying the 
theory in Refs. [T3] . [5|. 

It is interesting to compare the value EZI) with the analogous value, obtained earlier 
by the present authors (see p.393 of Ref. '§]) by extrapolating the nonrelativistic result 
concerning the positronium binding energy in a magnetic field to extreme relativistic 
region: 


Pi 


a 2 m 2 


h P cr I NONRELATIVISTIC 


exp 


2V2' 


= R n • 10 


164 


a 


(98) 


Such is the magnetic held that makes the binding energy of the lowest energy state 
equal to (-2m). (This is worth comparing with the magnetic held, estimated [2Zj 
as a 2 exp(2/a)Po, that makes the mass defect of the nonrelativistic hydrogen atom 
comparable with the electron rest mass. A more exact nonrelativistic value for this 
quantity may be found using the asymptotic consideration in [2B])- We see that 
the relativistically enhanced attraction has resulted in a drastically lower value of 
the hypercritical magnetic held. Note the difference in the character of the essential 
nonanalyticity with respect to the coupling constant: it is exp(7r v /¥/ A /«) in (1$7I) and 
exp(2-\/2/ a) in (1981) . Another effect of relativistic enhancement is that within the semi- 
relativistic treatment of the Bethe-Salpeter equation (Qj, as well as within the one using 
the Schrodinger equation jHj, only the lowest level could acquire unlimited negative 
energy with the growth of the magnetic held, whereas according to ESI) in our fully 
relativistic treatment all excited levels with n > 1 are subjected to the falling to the 
center and can reach in turn the point Pm = 0. 

Let us see now, how the result EH) is altered if the cut-off procedure of Ref. [l| is 
used. Consider equation EH) in the domain So < s < 00 , but replace it with another 
equation 


d 2 T 0 (s) 

ds 2 


4a | 1 

n 2 4 %(s) +m 2 T 0 (s) 
s o 


0 


(99) 
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in the domain 0 < s < Sq. The singular potential is replaced by a constant near the origin 
in (©)• Demand, in place of (THfll . that T 0 (0) = 0, ('h / 0 (s 0 )/T 0 (s 0 )) = ('h / (s 0 )/T(s 0 )). 
Then the result m will be modihed by the factor 


\f¥+l cot (v + i) - I 

which may be taken at the value a = 0. Thus, the result (EJ is only modihed by a 
factor of exp(—4/3) ~ 0.25. Generally, the estimate of the limiting magnetic held (19711 is 
practically nonsensitive to the way of cut-off, in other words to any solution of the initial 
equation inside the region 0 < s < So, where the magnetic held does not dominate over 
the mutual attraction force between the electron and positron. This fact takes place, 
because the term {^l' i ^ 2 ^/o^) ~ 65, singular in a, is prevailing in (ED, the details of the 
behavior of the wave function close to the origin s = 0 being not essential against its 
background. Although numerically the resulting value of the crucial magnetic held is 
affected but very little, we must bear in mind that the whole cut-off approach is not 
adequate, as argued in Section 2, and is burdened by the blind extension of the standard 
quantum mechanics to the situation, where electron and positron are brought together 
closer than 10 _14 m _1 ! The ultra-relativistic state = 0 arising within this approach 
is an ordinary bound state, not the confined state described above. 




4-2. Radiative corrections 


Mass radiative corrections should be taken into account by inserting the mass operator 
into the Dirac differential operators in the l.-li. sides of the Bethe-Salpeter equation © 
or ©. We shall estimate now, whether this may affect the above conclusions concerning 
the positronium mass compensation by the binding energy. It was believed that the 
radiative corrections to the electron mass are able themselves to annihilate the electron 
mass due to the interaction of anomalous magnetic moment with the external magnetic 
held. However, this moment is not constant, but becomes negative for sufficiently strong 
magnetic held |2?I]. So we are left with the primary result that in the strong magnetic 
held the mass of an electron in Landau ground state grows with the held B as [SU] 

fh = m (l + -y- In 2 . (101) 

V 4tt B 0 J 


For B ~ -^hpcr the corrected mass makes m = 3.45m. This implies that the mass 
annihilation due to the falling to the center requires a held somewhat larger than 
©. To determine its value, substitute m (USD for m and Lb = (e.B) l ^ 2 for s 0 into 
equation © with n — 1. The resulting equation for the hrst hypercritical magnetic 
held, modihed by the mass radiative corrections, B corr 


1 I 01 1 2 -®corr 

1 + — In 
4n B 0 


= 4- 


B r 


Bn 


exp 


has the numerical solution: P cor r = 13 P^pcr- 



( 102 ) 
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We state that this correction, increasing the first hypercritical value B}^ cr by a little 
more than one order of magnitude, is not essential bearing in mind the huge values (EZI) 
of the latter. 

4-3. Second hypercritical field. 

The same as in Section 2, we may attribute the 0(1.1)- symmetrical solution (PJT1) . 
which is a spinorial singlet, to the vacuum. Its 2-momentum quantum numbers P\\ 
are zero, which means that the vacuum state is constant with respect to the center-of- 
mass position x e + £ p of its constituents. On the contrary, with respect to the space-like 
interval s between the vacuum constituents the wave function decreases if these are taken 
apart, as stated above in Subsection 4.1; this means that these constituents are strongly 
localized. (The vacuum constituents may be thought of as delocalized in the ’’internal 
coordinate space” obtained by mapping the singular point of coincidence s = 0 to the 
negative infinity - see Refs. m, ia for a detailed explanation of associated matters). 

In this subsection we discuss in a qualitative way the situation that may take place 
when the magnetic held exceeds the first hypercritical value EH)- The eigenvalues of the 
Bethe-Salpeter equation El for the total 2-momentum components Po ,3 of the e + -e 
system are now expected to shift into the space-like region, whereas for B < B ^ cr the 
center-of-mass 2-momentum of the then real pair was, naturally, time-like. 

With Po ,3 A 0 equation E!1 becomes more complicated as compared to the case 
To, 3 = 0 considered above in this section. So, decomposition (1821) is no longer sufficient, 
but should be supplemented by an extra term Pfi f> 4 . The resulting set of equations for 
<h’s does not split now, unlike the set (1831) . EH) did. Nevertheless, at least for far space¬ 
like Tj|, Tjj 2 -C —4m 2 , the situation can be modelled by the same equation as El, but 
with the large negative quantity m 2 + P 2 /4 substituted for m 2 . Then the McDonald 
function EH) is replaced by the Hankel function containing two oscillating exponentials 
for large space-like intervals s 

exp|iP,5L±^}. (103) 

and two oscillating exponentials 

s ±2l V ¥ (104) 

for small ones. If one passes to the Lorentz frame, where P 0 = 0, P 3 7 ^ 0, and sets 
the time arguments in the two-time Bethe-Salpeter amplitude equal to one another: 
Xq = (Cq, one folds that the solution oscillates along the magnetic field with respect to 
the relative coordinate £3 — £3 (mutually free particles) and with respect to the c.m. 
coordinate £3 + £3 (vacuum lattice). 

We are now in the kinematical domain called sector IV, or deconhnement sector 
in Refs. uni, [si- In this sector the constituents are free at large intervals and near 
the singular point s = 0. The wave incoming from infinity is partially reflected, but 


exp < ±i.s 
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m 2 + 


p i 
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partially penetrates to the singular point, the probability of creation of the delocalized 
(free) states is determined by the barrier transmission coefficient [2|. Such states 
may exist if one succeeds to satisfy self-adjoint boundary conditions . These are, for 
instance, periodic conditions, to be imposed on the lower and upper boundaries of the 
normalization volume, in stead of the standing wave condition m, appropriate in sector 
III. The corresponding eigenfunctions are studied in [d- The possibility to obey them 
is provided again by the falling to the center. 

Now, the delocalized states in two-dimensional Minkowsky space correspond to 
electron and positron that circle along Larmour orbits with very small radii in the plane 
orthogonal to the magnetic field and simultaneously perform, when the interval between 
them is large, a free motion along the magnetic held. They possess magnetic moments 
and seem to be capable of screening the magnetic held. This provides the mechanism 
that prevents the classical magnetic held from being larger than the hypothetical value, 
second hypercritical held, for which the delocalization hrst appears. 

The tachyonic character of the vacuum state, i.e. the space-likeness of its 2- 
momentum quantum number, does not make a difficulty. The presence of this quantum 
number implies the break down of the invariance under the Lorentz transformation 
along the magnetic held due to the appearance of the lattice in the frame P 0 = 0 or 
of a superluminal wave in arbitrary frame. This is not a contradiction, since such a 
wave appears in response to a simultaneous increase of the constant magnetic held in 
the whole space, which is already a noncausal procedure. 

5. Conclusion 

In Section 3 we derived the fully relativistic two-dimensional form that the differential 
Bethe-Salpeter equation for the electron-positron system takes in the limit of infinite 
constant and homogeneous magnetic held imposed on the system. We studied the 
falling to the center phenomenon inherent in this equation basing on exactly relativistic 
treatment of the relative motion of the electron and positron. Thanks to this 
phenomenon, at a certain finite value of the magnetic held m called here the first 
hypercritical value the positronium level deepens so much that the rest energy of the 
system is completely compensated for by the mass defect. The most symmetrical 
solution of the Bethe-Salpeter equation corresponding to the center-of-mass momentum 
equal to zero may be attributed to the vacuum. In Subsection 4.1 we described 
the vacuum restructuring that takes place after the magnetic held exceeds the hrst 
hypercritical value in terms of formation of localized states of the pair, which are either 
’’confined” or tightly bound - depending on whether the theory of the falling to the 
center in Refs. [in!, [si is appealed to or not. We estimate in Subsection 4.2 the 
modification of the hrst hypercritical value of the magnetic held that may be introduced 
by the mass corrections to the Dirac held propagator in the strong magnetic held. In 
Subsection 4.3 we discuss the second hypercritical value of the magnetic held where a 
lattice appears in the vacuum and the latter becomes unstable under the delocalization 


Positronium collapse 


30 


of the states of the pair, the delocalized charged particles on the Larmour orbits being 
capable of screening the external held and thus setting a limit to its growth. 

The above limiting values are obtained within pure quantum electrodynamics. Up 
to now, it was accepted that the vacuum of this theory is stable with any magnetic held, 
contrary to electric held and contrary to non-Abelian gauge held theories like QCD. In 
spite of the huge values, expected to be present, perhaps, only in superconducting cosmic 
strings [31], the values obtained may be important as setting the limits of applicability 
of QED. 

As being due to the special, non-perturbational mechanism described above, the 
hypercritical held is determined by the inverse square root of the fine structure constant 
elevated to the exponent. This makes it hundred or so orders of magnitude smaller 
than other known typical values of the magnetic held that may be expected to 
lead us beyond the scope of coverage of QED owing to the lack of asymptotic freedom. 
For instance jj33,, the photon becomes a tachyon in the magnetic held of the order of 
Bq exp(37r/a). 
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